


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1991-09 


Performance of fast frequency-hopped 

noncoherent MFSK conventional and 
self-normalization receivers over Ricilan- and 
Rayleigh-faded channel with partial-band interference 


Lee, Kang Yeun 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/28447 
Copyright is reserved by the copyright owner 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
th D U DLEY research mate = and institutional publications c reated by she NPS community. 
«iit | : 3 Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
IN | KN Ox appointed — and published -- scholarly author. 
: ] LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


ee 

be af 7 yt) dake eye 
AEA +e Pil! cA iy 

EA Yell iat Ded Ri ve we Pa i s'y'g’ 
oe s : (yo able 
ae BN Atieah ied dy dit ', RAT eRe ARAN nae? 
p Alem Mate 962 
re! 


4 ofl ed ae ys 
4) ae ls Prem Petes fr. wnt oly vy a 


’ Nery Ws 
SHpyiot ot uaa *dyy ae 
daly’, aye Vial 

































































































































































































































































































































































































A LUA rea 4 Kas Vat a 
om 82 * oy ASX whe b +A 4 i 9. AD ond! hh Mie Hs £Y 1) tha AY ry, yes KH cae 
y ‘ he 9.44, Anis? } 1, *® s H ah i pe £ Re he! 3 25 a 55 ait a ety Fr’ : . yt 
: ' at UNS 44% Dee a) pag ‘ee ty Ao My: Ay enranees's Wh A Pete ea Ray’ gee Ga'y"ash 
re OMe | . 4 297" ’¢ A can A ghey’ vpy! i ie ue vy seo 9 rus A) o's a’ pa “aye ay 15 7 fo ® : ou ya 
4: $ , i sy ’ oe eg ef ae ¥ * Va ¥ | PrN A 5 "4 Wie ct Ia reatatt Beas ak 4, 4 lS ve sah * Me iy CAS i 4 v4 af 4 eee oa 
- ogy! i | asa), Ay ee "yf! » 4 M eat t fot 
yh! fy 4 ! 2 $ nat yt” Vee SON +9! if De Wh p oe sey tat adel dy 1, aah qi ‘ 2 RH As is ra 3 e} “see % 
et rt Ly | > 13 Poa a1 ya d wuP Ny, stat y vi ys ¥? TAN TENUM 6) yan mata 6 R) we SPO ie 
- a } Ss be ve y y's £30) tet sydde's aN ) dy’ ey eet tditasyl! + qa Ba i rte ‘ 
" hear 7,3" 3,527 HACE gat et ya! te nes : 
yt ne daze ettad yd ial, ni Peat BSS von ey cae Bie ¢ iP ahs 
Ave p et 1 U aie, bee Aw i sarsl Was ; ? he. a 
i yy yeas 13K) mt yf kart} x \% vey ‘ A? af ' ayes aay SAY ats sta A 3¢ fi en Me ' oe 
IAT By! SP St Ay sly yh? eels Hee SD RA De Once) vst: tye Ade ry EARS 
i i hO.3 j Pst oth Hs dra! ry vy ; ys MM $94 ash ready Ay a's’e iv », ee aay ik . er Ofte ae 
P: Os t Awe Us = yt, AOS }a® f . y " 
AS Oy te Daas dete tt Redes sil ehtiey enh atpanish Da Palakyty gab faa yee ey ore NN 
by8 a7 Tyrer’ vA? Ny ty Whee A) RA ae 9h, i i ‘ yan nee) Mn ‘ ae, mI Pet tae aie a , *, nae A 
“tee oy >¥ a", th iad ey h utd Wr, a A AAT rate as Salar (i mas Ae ab iptents a, jai Pic 4 a 
Tn Rath Aisa ae nae Taman Ma eet Shee KD igeoan a ee: 
Jp ey! bf Nee? A ee abet of ey’ ye Paes As pris rh ta 2 oy ah te yy ea hah ited Wi tem a me het ary 
A i) bee ' A “w F j if ‘ | , ‘ f y tye, th ¥ a 
, wir ty ty gtieg KA RAY ASS det) 5 hte 3 yf Wap d yy Waar pista, Os ie) nfote aes 
rt, gag Mp af ay + ay bgt gee Aa) of’ NICS Sh 4! Me, 
Hane hed ae, Qs’ Sait ih, Ke nha OSU nN vy ie oa Pai “i 
Pgs Ty) Daten aid eu pal de AL Eh a | ab Pata? “y"4 aang j i 
5 é) ee ete ‘, 5 nee b, ;: n saat, > | ® st 9’ yyy ys BAPE re ir ¢y aN wa es 
F ai Me pea om oy 3 Metab sie Paty IK Wy) I rs p % aca wate Sy AtStt § 5 ae 
? Fyty I j ‘ Pel a4) P ce 1 ye vy ’ cA ay _ Ke ie grat had Peta eZ i ] +4 a CERN vy B Mies fon | 
a ks j He wee | i GH . dy coda re CH Sata aa aia via 739 et 4 EN, 
ieee 2 25 Did ly, geateody hate y's Pik oh gen 2) Fy! at note fat pa? ya i. 1") Ne) a ely sates her 
5 » ‘yf: Ue i iy Files 5). ey 5 y “os Aye f. ShOTR - fie? , Wire Raby iran Baye! | Rh ae i68) wh olny pats 
x . 5 ; iP i Be Pe iat e set 4 ah a} Ata, Pabst ee ay sf gia’ s\a tah rp sate Vo Fi ada sinh mas nth me 
‘ ta aby iy hd RPA ne A ala Hag » pr A " bef Me r 
ie at v ay oes neat ey Ke Oh nih, vy" 41 pnt PN é ay 99) Beit 4 LAW, es te tee rD Cth ae ate abe ek aves 
Ps 48 i A te USA . + 7% aos ti Bag oFa'y * 
Pe RAR A cheAMy: tafser tsi viele’ inte) Sitter wie seas Isha! : ae ae 0% ao aoa % tise ay 2 
ytd hy ARaet ge lpnay yi ny ae MAS : Why , ath tet Tyla 7 Cetin Ri eas B us Lith Race ed oe gi : Rasen gives ih he Feat 
Bsa} ,f. he ue 9" ksh hig y 'yinde? fp Pyeng? or aten, + mt High 988. bed 74 . 
$ eetnty? “a | hs ya “y, ns iy siNry? Jy ahs Yah ty rah, 4 yf ’ i 9. ayer et ey SURO ae re} ae as a Ceipats| rer 
Oe Soe gud a eet shal aa “ae dria RaSh ie ghehity ly vegly aon v8 see Noe ey eae Se 
, me’ wf ; . tha* be } a r} ad me) 2 t ' i> s ‘ 3 j 
ys PERN A Mt 4 as, bas ot 30 ae As AS ay) wl ai feria aie A) tan as ei ‘iat ed me et } ie die pee my aS ng 
« 19's rr a i aay ist tat +, Ht pats te , paddy! Ay ” an Nese pi 4 7 if eas af y°%4. oe ; * ea) ie Bree “2 : ee ae ca 
Igiya %, ae ry ry ty fads ba nant? Wes My? a 3 edie pt 2b p H asf ‘ saat ae et oh See rst ® 
eas Pay’ ylna’y igtadaphtattety he er Ge at Weed Ai Setata’y a 99:4": , f 7 sie RSL Ne ¥ BEN ‘ nae Hae reas eI on anes tee 0 fab 
' % ‘) ; SRA “y ] isfah ats Wy Ket oY La ‘ aly eee ry* i a ' o ane # ahs aeae 641 ekcee ri i o at Pa se re 2 
f Pa Why ya 1 rh Agila AUS phat ), ‘ duty? Bead af hn vi) he % + ane Vitae nt ‘ee Ae "7 aes sit ne {3° sla can MY, AAs ao Hob o 
doe: Seren He) 46 HEP Robe a’ Wee n i ny byte Obs ie gn Aa hParge of ’ hei at Lat ies 
Pryde v2 PEE ROE tt cay Bo - alate CE AREF Rae ee x ou vs one spite yapasahe 
MAT ar hye Bal  aly thet Jaragial y af. Peadaty bys ate st 9) Cree sad Lf was tke Y ‘ aN . 4: 
7, yaly’yy's iy “¥ yo} ay hintal ye Ren ig Sante pil ys id jou citys ‘ Vijat Fars : Soe br he wary 
bd, Da hat oF ty aS MIRNA MR A RMS REM ON CR 7 é ees i Be shy 
‘h piyis Ae eal us M natin e VP dtatatghy thd WV. } Ay + 4, * et Peent » righty a one ' fos eet . Rng 2 Sf ca WALA Tatar " ae fa? ad ee 
yatety hale Betty lye?) yey LS AAACN § np wy cose Ag AAAS: i Rati poy tet doFe Sane at oe 7. ; Bsa br 
2 yy, A iB Ag vy yiat “it 3,0 a hs. sah! hs P HY AK 1,3 ‘ Ae Ee tar 5 ' ON ys vi ue hee oe a " b sak rh 
os ae ote EP naby gigladytttys her ae sesh qrpls? 4's 8 We EMC ee ~ 
5 vy 1) 8 ee any GS My yp Nee hee My % ,; Hit wrel 1} ; g }. r 1 te b2 a Hy Aa A sare iar) ie ig 4 
ye bind ae IID aad ob yeyiaty zy? eye! dy Te, Ay ips ata 4303, ney * Patines An Pr Rana 
FLUNG COREA RS he teh eet Hohe di fab plyly ‘3 Ma Fi dgtdi ty beta alg IMCS Seay aire ak ae Bit Lis 
2 : sits bist cal Th Sf Lae [hat » y +h, ¥, wat ; Te sis A jiae yyy ta sry : fi Pint "a3 He, oft fe Py RES art seta as by 
te eye ees La le aie pela ear aay sts edt) Ayia ei Lee 
a 9 e $ ' ¢ J 13 ? e ’ a te ys 6 4 tus hie 
‘et yty Avie i i Ba: A! os Se 73? oe bat Ay eat? zy er: 
: aly Y Pa, fi > ; Af 4 ed es 
ae : Wie ft hy? at als : ee re 
'¢ > hy Tits i. t 
f } a? # fi i ‘yk vipa yi at ei aleyas ; 
‘ ‘ eee y es, ve & Pad vs 
: + kK 1 “3 I Ars, Phe Ty x at ie sf 
or “fy Fy RM ten. mats rit yt 37 *, a be ea ? 
‘ 1% 8 rh : 4 ; & Aoi, eh 9 F Ae iat ty 
Le eS ee ae ya? if Meat vpatyl 72 a jit Ria 4 ey ay sith dats Rey 30 
SEINE RAH IR eae nae vy a NAMIE, a he ae OSNRee aaa a ae i 
iat, af BIE ye 32 Aye ees : Rg Pda t aid rshohs , ys! f shee want rast ake Oty oh Mytatnd me neg? i Sey iho 
Me ee atu sab datntet oth ee pale at seecary yl at Nae . ity AV. ys Ny ?ahy fast 
m # Pe A . hes ea ee Poel reas , eA MO wR % Tad ” 4 oA: ty! ie < yes hte i .: 
, lh Be fee eri ye URBAN Uteeias Pydy, teatiy ‘ an a 
2 Get kA at RSS I ln Per SHEN MR oe Ben oe 
vm § Se pian) tat abel chy wos CRP a aati Be Pe ae ‘ x«: tab Pyt nga dalal aC ath ty gt 
4112 A) + HD Va had’ ‘ RRC OR ARC nie rye } GS : *% ») niguiite nat seen Kt: 
eS ee Wa Ya RL Be aah FMI Gehan rate wnt he 4 ear ia H Pee a Sail ist 
oy Ts 3 y,/, Bh, 4 of. aA in: an id out AS, fais! vy? WE yy 1 wh ‘ef iY tas ete ' Ay 3 tn? a st ee ie it ade A: 4 v LOA: 
I / p bi rH Ae Sais: eh ’ ay vy Cia sy _ i 14, tar iiss ‘ ; OI % 
tt, y is rs Pee U3ix a Lb oe Ul b ¥3: sty ee PRAY fy fy ‘ : i ayte Pad obi Peale 5f¥ vay Pera? Ata 
hid a tate ye ‘ "ly Py DA A De fat igi ve! ashen nis Ne es ah i A ere By “need YO” in hate a oy? oe : Ly eae fiaMe 
Phe ae aT ee ee ey ATE Ny hips 4 et SNM TMU ed RA ORLA ae AEM Re Soe AEN us ae Neon. 
4 Tate “Th ist a ray Roae Sh hd es oe 2 ak : 9 A, iat Ae Raa f ‘ ; ies Matar v nt a ile? ata ¥ ‘ Ait S083 eye ey oe pated ae 
Ng hl ie FAO of fy + ie es Fad y's we cule rs d aS Bis . $ he ¢ Be o™ SP yk - ta be) ., ed syinen 
~ . snes Por sbet tated Tyler ty sate Me i ‘ rae “i vy wes ar \ WN Ben 
Re yabstes j a ee ay te wits bad Mig iy , ' Aig! Fa gta) ta Augvsuent RODE Fu racvat Fat ideas a 
Aa Pah « Asse . eo ee Rees na HIS : ak ; WeNWer eda ; oy yoy tate Ba ee ee sie 
EEA ey OTR orc Padshy Bel Ae Satta ary daly Lhe ° Sete of F Pe! ‘ Se Faye 
Pate ger 4 tyres? OE as thie aah mee Pyle yyy rygns taba tes a mae ae) RAS a 
Paar, be eT une bByd fyenr at yh; i ' "3 Q Fal ry" " us ky 
Fa pPrxin a, 
dont 
. 


a2 Bee na G 
ag oO Gah Felts Eade 
iecane, 


Ag? a8 
oe ms oP: ane pitt 

































































































































































Byatieas eS ee 
C0 > Fe : Map iy ee eit st 
ty ATE EY be ry ee oe ee cory ite See 
van De a BN ee oh gh HAL Hy haires! een yo ay, ee 
2) .. Fob pe ee LATE 8 ; : he ay yt 84 ae ; ‘ er , rev F Sad een aa ENA PCL See este Pred pees 
. : y : y} yy ie iaaaaie Ya Cyt 7 He Sap GER bP MY . 44 be id Chee peak watly iets d eee VET Ses rt Sees 
ED lel AR Mey ke SPATE ML TAME Wt 48 ed e, a es ‘ Pays tee eae cts Be Bece Msn Pet PONTE 
WP ta SiR Nouba bake TURES y Ma eS HeaR SRNL ME SMES ORAM eee nana tae ae De Onin e nee aun ceen ey ae 
b eee Cbd tte Fry Rae ’ ii; 4 | ‘a@ > + rind fe 4 PeRTS \ Ke 4 a A ? bps +f * 
. hs ae i ! ‘iy st as ae sty ? stati dybyh 2. Aire vy \ 4 ye at ry Mrs df 4 4 Py , 4 5 ; 5 me 7S} mee ayy pee ve fy v7 . ant of ts iz UD eae eas 
- ; : jis Ay +4 ' « ay Ghats” ra ide ‘ty! B f H ‘ b-¢ ? 
ehised tas ou yess aby 1 tina MI KD hay ys Me Ries As yey Pel BALE gr aligrah yon ? oe Ta gy fe yea W865 Pye tite rat ve Re nis ef 
yu yen Nt5 Bae yraltey pe? sone BASS nace SE ; AH eas 2 rd fey ne ie) ta 3 ie ce RS Tat hee een telat 
Jot hy Ae) Ma PRA? ST etpraky lat ye catel sta ty et ‘ “ ee AS tye ae eS ay Tah a eee Ly ae Tey ieee | ae af LEA CEN 
9 A ? Hd fie a's 3") ' 84 2 Wining Ieee % iis vj 3/36 gis7s!s ray Oa las Sig aay A es eb ie RAR 
ae aa at Cae ORES RR POS Te hs “ Thy. 
ie woke ie Wivae eee eit ey brs , ad vent 
! HP OG, hy ore % ) jet. i) Ryo ; ifn pyre AF hint i: Loti ia rN at wie) ? § Et is Aged en we sth ne x 
ibs o8 tos RAS t ry yy a sa EA) ‘4 an rey Sins Cy Py) 0? asi 5 5 +f, oe - eh 
waa yihFe q ys} : Pa? p 
} 173 3 iy Be aA PF me Tae FH ey: tv ihn o, My! Als ty f oe if “ . f oe a! ne hee 
i Oey } ers ;' eat er ry , ; ahy bie U1 radar ‘abate a' Nai , Ge ) ‘ 
ry abd, F we fi Surya 4 yt eee Py poy ** iat F) at, Anis 3) Nhe : ir Age 
ea bed, t Pisa, Jy 179.3. 8 Fda ra = i Oh Het Paty f ists 
vag of dye a a bi a it et Sa Ine lb aide ay po}? oth a oe ies ete aya iaes 
roar, FPG a, les _s iy a} MW idy rei 15! iw. f yf : Pre ‘ , aoe P S rt ? ses a4 +) ee 
a fae * FG. } ater | i Anan! Hil yas ray A ibs ; ‘4 54; Th * ayhy rst . Pa e nee ad pee 4 <7 ate Hie 
pene yt py pe bhi eee ee se pe pe For ey Ars Sethe fi Py #yInt bao a fs 4ye Poa ¥. 
baniads LAMM hy | i, Set aN ep) Ret CLE Ate iy Mn Tice thea aye tan ty. ita ee Ae Rete 
t ah a byt Pon i ped 4 BS eyed ur Nes as sf : ei EN ix} apa Pit a ‘ fire a Site Ea i AE y sry ee ot: 
- a $ 4 Ase at ’ § at FA yrs Ete sy thy fas Lous A ate . " J ta ¥ \ 
Py ¥ dys 3.) ae y's "5 z Me vie inf i ‘ci Aye i, SO ae fy MIS pouag Naot Bay f RAS 5 ry x a, : toa e iriets a 
a p y } sey i Mae mehad Poe Bd cad fae ¥S bie Oe MIEN Egan he : nas * Mid Biss ishe siete gt oe #: f st ot , 7 & ees, 
Wu Hah kh feaaght : Pade n} Ws o Lat tails ey SRA Wale i 3 i sees) es | me od ‘ * rf #8 f : . st 
ie map. Ge ! ‘ VI yh ~ fg ls Pad. nt hm dads oe al . hs 
i yre yer | he yy pis)! aad ea ie My a." » iiyisobs hag tay satyh ' ‘ 2 f.? x 
bP ape tae a's Re Ye «fay sfeeetaand Si 43 ; 
; fae 4 “Ast "ate pa i gis) as Wht ae i ite rhe sty 12? ry re 
' - 4h 8 a) ey +02 Sate nts a hy . . } ea ; 4 MF 7 
‘yy. ved a! Pb |) Nasi gs ‘ As yan Ahi i ; i¢ 5 ai Ane (y i 5 : Pat. 844 3 fa! a Gis 
: Fad U4 e* syd te 4 Beu ek : ee a is i. 4 Ls Pe Rit ty rai] Reds iy : bat ON : Bae Ay ( Syie ea wf 
i i ‘sma vat oer ee | } pete aty 2 ‘ ‘1 aed pe v AY. Ned SNS BP ee 4 A, ie if dhe BAF toy? al yc Bas sof Ose xt PF Rg . ’ 3 . rage 
par ie: Jas ; Vt v ae a and Aa AY 3) nye) 4, ieee ish We we v ‘ uae ¥ ot ‘ pete fi ne i . ry $y f Sy 
uae ar (no ae | TS SAE gh ee cB A et edits Hehidy te 20 2 : . at AGS ; Mei had ati 
aL iri * 6} re PA eat, yh gh Re sais} FR pls nse AS yp ed bie Tat * Fag ied: ° 
My a pte wats Fo gees ti Ween rihh Hy si I TAME " ASAT NON EH ff Pays: WKS y at wae j .* ei ’ 2s 2 eee Y , 
7 ae i% ON wan? ? ZL Sirlicey 1: 4 : ry a sfite rt i eH) wes Cae] oy ak yi Lae i ¢ " 
) id re 7 A i as auw it j She | i ty, PA Pe Pe he ee} yt ry tus Ahead hee = Sate 
lefauae a il Bercy, Eh keke 4 je Pd iy ta, Sule Pi 4.8 Ae 2. r4y pF is feat Anes ty , Y ates 
Te Dra A CHRO take yah hh), Bb atmd eat? 4 PRS MAE, peat, Mahe eity fg ar oala a! 
sgl 7, a) GU tea or et ee ee what be ae oe ty! ae Yu Hee ' 
PUI We ik ad OP De oh O08 ale, CaM 8! Paging wt) A ee ah At) 4 ed | ‘ eye Reha tere 
, jy? a 4) Nartege a Ft \ et Btn ris '’ hae ‘A Ny ele au ae BH 7 ‘ { : x 
ri hy Y nel dD ge ‘, ; high i) ¥.€. ie) ry Pvt tee) Ae Seah Ha TY Mee h dig. eat } e Ay. Fe 
hee s* LS | A0"} Mg — ei } ee ed 4 He IF} pra ¥ ee te 4 a 9 Tear 
Me ba her ho yy yVvG he Lee oid we ih eae Fe tavye gly ip is H a3 H <4. ef 
4 gyi es 23} 4 fs ; x he ta PORTE a hs eS RS k hte : es * TP , { ' ray + : REM ie ¥ 
f ry} vuln ii f H BAT, ibaa 545 Sova ihe a ily; Beets a FR. . ty ? se eye reid oe ut 
. “Y htel ets he a ; t yy nt uy PAW mae ee mip ited vis ie a) } y AE) ay Iie 
' y 4 Ri: 4 , ihe \ ~ ; xa 
; 4 bs are Ne 5° al iM 4 Ar ye Fars ° 456 J 's 5 Ff ‘ iy AE A) Pe , Nha dyytdes EAeNs ma rege 
uf ; Ji ‘. t y af er ’ i$ i s ate, ht, A, y ¥¥ Aaw) By <a ae ¥ # dae, 
, sy ot Fe yeas Wad ne nae i. Atte as Si dasis whey, 
CAME MG RSS he MESH aT te gO a haha Iai RNS) 
a ’ “i MWe AE He , ‘sf 1g <i> ¢ yest hyn 4 Sen eee a eye v4 9) 
ea br wa Gbee yd KOS eae Paine: tLe ital, oa, 3, 26 Sus te aay ey. 
gh i? 2G " i aed fis io We, HDs Eade jt ii YY oy ne A‘ yet Aly ashy WY “i 
3 ee Fee? 1 a ay tye | Ae, Td hy J ara ¥3 4 hier */ 
‘ ’ ‘ss 4 ry 2° 4 , hd, mre? iy | i ' A us vy ins. j 
“a! ht: j Ae La fo ae *4P)4/ a 43.3 3 £8 tet a an 
hi’ : , sh eT e. 6oe ie ’ 2 4} Bigs ay AL Le Re A 
i 7 Had rr he oe he tit by ye © 6% oH {et 1,9 ee f AY Frond. & FY to 
17 Ae Ee yi Gat eae Bah Met Nhat DET Tae Rte facut a oH 
MET 4 Yet.¥i oy yt a Vid Bia he a hry hess ‘yf att e at 7 Ray Ata) x ae ay ay te nein " Rea 
us i) raya 1 2Gh fo ¥ Ae Re ee bw kgs vee Risdon mf a, vp hang. pes a ae in Lesions Hig * rt ¥. Hes ey Ce ® ride 
* yk saa ONO ah aes Bh Wea ik fe i Oe cael Se oh Wee 
os Fy £43 Oe Bib & a ie) ny x az AAS wh . ue 4 ‘ i q i voip 
; ew dss , 5 4) yf ¥)21, ee dy Ty . iS, ifs ; i ZY ty eA wi he, *«. iy TAS 3 Hy: Bie host sf ¥RS aa 
7 ; f Ve Nt he q bint AA 3 “ys : chy APs Res HOY, v8 \ pana 
é pips Foye Los bps * fi iN! Sh ad, ap ‘ ; ik: 4 J ie an Oe _ i sate KS ¥. 
. ae : ‘ ‘ 4 ‘ 1 ae Ler Poy foe a sey 
| aie NAPE tee FAD, ph Syybsety Dey aries mise we BASS ERE beth 
D ,¥9 aoa : da ’ rat th fa TH $A 3'; i; at 4 oe Seine Fre , 
inf i pe a er . 4 soe LP A, , ee eu ee i yy Pee ¥ eit é 
au ‘ eu "5 a Uy ' ! i be >> 5 + i] - va }. ' ty } & 3 3 es 5's Oi 
a ite le heat gia Nandita he. FE 9 4 
LPP ter Kit a a rult gh ieee eee bra” Saye 4 
: Bytes yy oy tg t Cry» FGA i eye Ts ‘ ? 
mt Da ae om Ark Poe PF Fm HE OL tat teeate aiylalg 
al? aid 5 mat > 4 gVin a” ot ao ayy ye ry j i > 
4 i irs i peg ie Ae acy Je, , i. ‘- 4 ys 
u% ee raee Be, Mat tir te t Meath vie ip ale at ai ic 
,,! \ iAH? ASI YK) Mt a ee abit i ] “anh a2) aeites Aa 
“i ! 044 os PY a a4) 7 ‘aed Yy ie Ea 3 , ar vant i Aaah 
Pele Be . e ha - #0)? Fy ' ‘ Ware a bk ) Aye at H TORK s/ 1 
. r ‘. ‘ ’ tet of P id! . yy F re ed P 433 4 4 Ki “9 hy ane hs xy wi af 
“a ‘; . ‘May Ay FOF 8, (? Th. } ee é 4 "y y Ret i 5 a 
Ln i | 4 ‘ rane . 4 a * 
ad : f ; ' Yon,” ha" ’ u% 4 : % J ty rT) Ce ¢p ; ae re ite oH iy 
KANON PON CAN BS PASO Dae iso ts Ante vat +r; hb? 
% a! bie 5 ie) . ? ae of ie * .? ‘ny xy « A, ) ag i . re a 7 24 , rx ES " 
not HI Re SER ON Phy AALS SON, hls toe te aM Head) we el. CAT ee 
ray ") are, 44 pa ee ey va a) af ye fees y's ‘ U a) UA ey ‘. evs ‘ Ue { ee le ies Vo 2) 2 ¥- 44,0) ) a » 
SUE "rs? LU 9h Thy 4 Ca te a » ) nt iy \ 1,0 £0 4a, “4 x r} - 5 st at ug OE ao Pte’ st 7 =* : 3h at ae TOA ye Fe . y 
es Or ht: P fe’y A}: Pe 8 > " mee ra alees CNSR NE pe + : Sate on OAS pea ate e 
x i 1) PPM ANE dw hy SY abi “los AP ee APT Aa OS er hae ee 5 1Ne ras rh Ne AN: 
i y aie i t Phe FM ty GF 44), og sh fo | sh bt 25 Y bi PoC eee y hy a AK 4% 2)" . a Ae | FS 
}' Vay 28 rT Ch. AS eae) a a ; Yi 7,9 PB mee } RAZ Pag *a°x a yt ae Se ate & F9.9é Rays 
‘ heats fos yay ian) af , +4 she Gea ON yea I \ y Wat es, a Aye ihe &} Ate FORO 4 Ey rks ‘ a 
8 ayy ‘ Pe ee a | 1. i ww, a) th REY Fe %5,, & MN a sR? ge > ak Ed hs Pi 7" @e i AY 43> 
UNL Rie ON MRC ae CMR ry A TCS hasnt a, RA ABTA VOR x Bie 4 Be HS in Masy 
NURS WE hh he “es AT ae Reet MCr pats g seh 5 ek RAT vf) atest Ae al AN SR tleay) 
a as oy hgh it ‘ 2 Le h Ns D tad ; "P . i! $39 . 3 San 3 Ae KIS ain! 4y NS , < ne gays ‘ is se) Le eee ae 
H ’ H_F oe i We . % Lf ' t ” | ws : bet S ua 
Nye He ee ah ier? rh 3 ’ rs ' yah. Ay rss aa fahs ait 2, J Sf ob), * —~ §- , hy fe Ce a 
. Ae ay i) i nA? ff 3 e 4 oe hn Sh ae Les - od sk % Ky ra. a ay thes nes LAL rh Phe 
¥ re Ol h PUI Ss U ‘ nk foe i. Na er YY 55,518.69 4 tos Shee ie axe i! ‘Fy A Py qx ’ete Ae Wa NA, 
’ an i2  nT a ee beh ERK ats: > Y ae Oe at. @ i et) Fe he id 2.009 rte WAY Be . 
ae a as t ‘uf ry) ‘) ry, i ? ty ae of Rs Wea "Oe SNe ae ASE a ate Sig ah wt REM q ‘ 
“oY Is, Tees a TRY eg et EY ie PO NN ee ct ra Lots « 
: hie ae ; rs roe | Le if Pad Pi bh ae ft teat sy AAS ‘ ie: 45% yal Ad Wi 33,4 a? Ap tht eA) | 
fee 2 fy yey J ss. 4 ye 11 Aa? ela FaVay se a Fo eM er ae a ae MBS ne > ri | 
pane, al ,',’ Vek EM pp ey £ eRe ahh! vi A fdas ge Ride. 
9 F445 ' : A ae nf t& 3 ; r) 4 
Ror ahhh oe See atta han At tat 
y iy ee Pe a | 5 04 cc 














NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


Performance of Fast Frequency-Hopped Noncoherent MFSK 
Conventional and Self-Normalization Receiver over Rician- and 
Rayleigh-faded channel with Partial-Band Interference 


by 


Lee, Kang Yeun 


Thesis Advisor: R. Clark Robertson 


Co-Advisor: i og Cad eee = 





Approved for public release; distribution is unlimited 


1257832 





classified 


urity classification of this page 
REPORT DOCUMENTATION PAGE 


Report Security Classification Unelassified 1b Restrictive Markings 








Security Classification Authority 3 Distribution Availability of Report 
| Declassification Downgrading Schedule Approved for public release; distribution is unlimited. 
Performing Orgamization Report Number(s) 5 Monitoring Organization Report Number(s) 
Name of Performing Organization 6b Office Symbol 7a Name of Monitoring Organization 
aval Postgraduate School (if applicable) 3A Naval Postgraduate Schoo! 
Address (city, state, and ZIP code) 7b Address (city. State, and ZIP code) 
emterey. | ce 939-43- aon Monterey, CA 939.43- 0) 
(if ipicctiek 


Address (city, state, and ZIP code) 10 Source of }10 Source of Funding Numbers sssi‘i‘itsésés*™r Numbers 


Title (include security classification) PERFORMANCE OF FAST FREQUENCY-HOPPED NONCOHERENT MFSK 
MmNVENTIONAL AND SELF-NORMALIZATION RECEIVERS OVER RICIAN- AND RAYLEIGH-FADED 
MANNEL WITH PARTIAL-BAND INTERFERENCE 


Personal Author(s) Lee, Kang Yeun 


a Type of Report 13b Time Covered 14. Date of Report (year, month, day) 15 Page Count 
faster’s Thesis From To September 199] 90 


Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official policy or po- 
lion of the Department of Defense or the U.S. Government. 


Cosati Codes 18 Subject Terms (continue on reverse if necessary and identify by block number) 


eld frequency hopping, spread spectrum, partial-band interference 


Abstract (continue on reverse if necessary and identify by block number } 

An error probability analysis is perfonned for an 1f-ARY orthogonal frequency shift keying (MP'SK) receiver employing 
st frequency-hopped (FH) spread spectrum waveforns transmitted over a fading channel with partial-band interference. 
he partial-band interference is modeled as a Gaussian process. Wideband thermal noise is also included in the analysis. 
versity is performed using multiple hops per data bit. Each diversity reception is assuined to fade independently according 
ya nonsclective Rician process. A nonlinear combination procedure referred to as self-normalization combining is employed 
y the receiver to minimize partial-band interference effects. 


Work Unit Accession No 


















) Distribution Availability of Abstract 
J unclassified unlimited (J same as report CJ DTIC users 


21 Abstract Security Classificauion 
Unclassilied 







1a Name of Responsible Individual 
. Clark Robertson 


22b Tclephone (include Area code) 22c Office Symbol] 


(408) 646-2383 ECR 


BePFORN 1473.84 MAR 83 APR edition may be used unul exhausted security classification of this page 
All other edittons are obsolete 








Unclassified 
I 


Approved for public release; distribution is unlimited 


Performance of Fast Frequency-Hopped Noncoherent MFSK 
Conventional and Self-Normalization Receivers over Rician- and 
Rayleigh-faded channel with Partial-Band Interference 


by 
Lee, Kang Yeun 
BSEE(Korean Airforce Academy) ., 


Major, Korean Airforce, Korea. 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN SYSTEMS ENGINEERING 
(Electronic Warfare) 


from the 
NAVAL POSTGRADUATE SCHOOL 


September, 1991 


ABSTRACT 


An error probability analysis is performed for an M-ary orthogonal frequency 
shift keying (MFSK) receiver employing fast frequency-hopped (FH) spread spectrum 
waveforms transmitted over a fading channel with partial-band interference. The 
partial-band interference is modeled as a Gaussian process. Wideband thermal noise is 
also included in the analysis. Diversity is performed using multiple hops per data bit. 
Each diversity reception is assumed to fade independently according to a nonselective 
Rician process. A nonlinear combination procedure referred to as self-normalization 


combining is employed by the receiver to minimize partial-band interference effects. 
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I. INTRODUCTION 


Spread spectrum communication systems are designed to use signals having a 
much wider bandwidth than ordinary communication systems. Some of advantages 
of spread spectrum signals are the energy gain achivable against a narrow band jam- 
mer, low probability of detection (LPD), low probability of intercept (LPI), multiple 
access operation, and the capability to overcome the effects of multipath [Ref. 1]. 
The primary disadvantages of spread spectrum signals are the time and frequency 
synchronization required by the receiver which make it difficult to implement the 
system and the wide frequency band requirement. There are three primary spread 
spectrum communication methods: time hopping (TH), frequency hopping (FH) and 
direct spreading (DS). This thesis presents an error probability analysis of two fast 
frequency-hopped M-ary orthogonal frequency-shift keyed (FFH/MFSK) systems. 
where the term ‘fast’ imphes that the frequency-hopping rate is greater than the 
symbol rate. The first to be examined is a conventional FFH/MFSK system and, 
the second is a self-normalization FFH/MFSK system. Both systems use noncoher- 
ent detection. The channel is assumed to experience both fading and partial-band 
interference. The performance of the two receivers is compared with each other as 
well as with another type of FFH/MFSK receiver referred to as a noise-normalization 


receiver {Ref. 2]. 


A. Conventional Receiver and Self-Normalization Receiver 
In the conventional FFH/MFSK receiver, illustrated in Fig. 1.1, the input 
signal passes through a wideband radio frequency (RF) filter and is mixed with the 


frequencies which are created pseudo-randomly in the frequency synthesizer to dehop 
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Figure 1.1: Block Diagram of Conventional FFH/MFSK noncoherent Re- 


ceiver 


and convert the signal frequency down to intermediate frequency (IF). The bandwidth 
of the RF filter (W’) must be wide enough to pass all possible hopping frequencies, 
while the bandwidth of the IF filter (B) needs pass only one hopping frequency. 
Thus, the bandwidth of the RF filter (W) is approximately equal to NB, where N is 
number of the hopping frequencies. The receiver has M different branches for each of 
the M different symbols. The IF signal is demodulated with a bandpass filter and a 
quadratic detector. The output signal from each quadratic detector is sampled L times 
and combined to provide a decision statistic between the M available outputs. The 
transmitter uses fast frequency-hopping to provide L independent samples through 
frequency diversity. The receiver is assumed to dehop the signal with perfect timing. 

In the self-normalization receiver, illustrated in Fig. 1.2, each of the M quadratic 
detector outputs are divided by the sum of the quadratic detector outputs from all 
M channels. This normalization factor is the only difference from the conventional! 
receiver. When partial-band interference present, the hops that are corrupted by the 
partial-band noise jamming dominate the performance of the conventional detector 
since each hop is equally weighted in the overall decision statistic. The normalization 
factor minimizes this effect in the self-normalization receiver. As with the conven- 
tional FFH/MFSK receiver, the L samples at the output of the MJ quadratic detectors 
are combined to provide a decision statistic between the M available quadratic detec- 


tor outputs, and the receiver is assumed to dehop the signal with perfect timing. 


B. Fading 

For both receivers investigated in this thesis, the received signal is assumed 
to be a combination of a single, nonfaded (direct) component and many reflected 
(diffuse) signals. By modelling the fading as Rician, in which the signal is considered 


to consist of a direct signal component and a diffuse signal component, we obtain a 
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Figure 1.2: Block Diagram of Self-Normalization FFH/MFSK noncoherent 


Receiver 


general result that is valid in the limit of large direct-to-diffuse signal power ratios 
for Gaussian channels and in the limit of small direct-to-diffuse signal power ratios 
for Rayleigh fading channels as well as the general case where the effects of both the 
direct and diffuse components of the signal must be included in the analysis. For a 


Rician faded signal with amplitude /2a,, the probability density function is [Ref. 3] 


ak a? +a? apa 
ae eae (sett) lo ( a u(dx) (tgs) 





20? 


where a? is the power of the direct signal component, 20? is the power of the diffuse 
signal component, u(e) is the unit step function, and Jo(¢) is the zeroth-order modified 
Bessel function of the first kind. When the direct signal component power is zero, 
which occurs in times of deep channel fading on the channel such as when the direct 
communication path is blocked by terrain or other obstacles, the probability density 


function reverts to the Rayleigh probability density function [Ref. 3] 


ay, az 
fa, (ax) = ag! (-3) u(ag) (Ale) 

Also in this thesis, we assume that the smallest spacing between frequency hop 
slots is larger than the coherence bandwidth of the channel [Ref. 4]. As a result, each 
hop of a symbol experiences independent fading. In addition, the channel for each 
hop of a symbol is modeled as a frequency-nonselective, slowly fading Rician process. 
This imphes that the signal bandwidth is much smaller than the coherence bandwidth 
of the channel and that the hop duration is much smaller than the coherence time of 
the channel [Ref. 5]. The latter assumption is equivalent to requiring the hop rate 
to be large compared to the Doppler spread of the channel. As a result, the signal 
amplitude can be modeled as a Rician random variable that remains fixed at least 


for the duration of a single hop. 


C. Partial-Band Interference 

One of the applications of frequency-hopped spread spectrum signals is to reduce 
receiver performance degradation due to narrowband interference (either intentional 
or otherwise). This is accomplished by transmitting each symbol at a different carrier 
frequency according to some apparently random pattern known only by the trans- 
mitter and the receiver. The bandwidth of the frequency-hopped signal is generally 
much larger than the signal bandwidth in the absence of frequency-hopping. Conse- 
quently, in order to interfere with the entire spread spectrum band, given a fixed total 
interference power the noise power spectral density of the narrowband interference is 
reduced. Rather than interfere with the entire spread-spectrum bandwidth simultane- 
ously with the consequent reduction in noise power spectral density, the narrowband 
interference may affect only a portion of the total spread spectrum bandwidth at 
any one time, and the portion of the total spread spectrum bandwidth experiencing 
narrowband interference may also change in an apparently random manner. This 
is referred to as partial-band interference. In the case of partial-band interference, 
only some of the transmitted symbols will be affected by narrowband interference, 
and the question arises as to whether performance can be improved by implementing 
diversity in the form of fast frequency-hopping; that is, each information symbol is 
transmitted multiple times, and each transmission is at a different apparently random 
carrier frequency within the overall spread spectrum bandwidth. 

As mentioned previously, in this thesis the performance of two different fast 
frequency-hopped Af-ary orthogonal frequency-shift keyed (FFH/MFSK) noncoher- 
ent receivers is investigated. The performance of the conventional FFH/MFSK non- 
coherent receiver in a Rician fading channel with only wideband noise (no partial- 
band interference) has been previously investigated [Ref. 6]. For the self-normalized 


FFH/MEFSK receiver, system performance is obtained as an upper bound on the 


probability of bit error, while for the conventional FFH/MFSK receiver, the actual 
probability of bit error is obtained. For both systems, the transmitter 1s assumed to 
inplement 

In addition to wideband Gaussian noise, the channel is assumed to be affected 
by narrowband noise in the form of partial-band interference. The interference that 
is considered in this thesis may be due to a partial-band jammer as well as other 
unintended narrowband interferences. The partial-band interference is modeled as 
additive Gaussian noise and is assumed to corrupt only a fraction 7, where 1 > y > 0, 
of the entire spread spectrum bandwidth at any one time. This is illustrated in Fig. 
1.3 where S represents the signal power level, N;/2 the jamming noise plus thermal 
noise power spectral density, No/2 the thermal noise power spectral density, B the 
bandwidth of one frequency cell, and W the entire frequency band. In addition, 
partial-band interference, when present, is assumed to be present in each branch of the 
MFSKk demodulator, and the fraction of the spread spectrum bandwidth experiencing 
partial-band interference is assumed to be the same for all hops of a symbol. Hence, 
y is the probability that partial-band interference is present in all Af channels of the 
receiver, and 1 — ¥ is the probability that partial-band interference is not present in 
all Af channels of the receiver. If N;/2 is defined to be the average power spectral 
density of narrowband interference over the entire spread spectrum bandwidth, then 
y~'N;/2 is the power spectral density of partial-band interference when it is present. 
Thermal noise and other wideband noise that corrupt the channel are modeled as 
additive white Gaussian noise, and the power spectral density of this wideband noise is 
defined as No/2. Thus. the power spectral density of the total noise is y~*N;/2+.No/2 
when partial-band interference is present and No/2 otherwise. If the equivalent noise 
bandwidth of each bandpass filter in both the conventional and the self-normalized 


FFH/MESK demodulators is B Hz, then for each hop the signal is received with noise 
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Figure 1.3: Thermal Noise and Partial-Band Jamming Noise Model 


of power NoB with probability 1—7y when partial-band interference 1s not present and 
with noise of power (77! N; + No)B with probability 7 when partial-band interference 
is present. 

The symbol rate is R, = R,/ log, M where M is the order of the MFSW modula- 
tion and /%, is the bit rate. Since the MFSK signal is assumed to perform L hops per 
symbol, the hop rate is Ry = LR,. The equivalent noise bandwidth of the bandpass 
filters in each of the Af channels of both the self-normalized FFH/MFSK receiver and 
the conventional FFH/MFSK receiver must be at least as wide as the hop rate, and 
in this thesis we use B = R,. The overall system bandwidth is assumed to be very 


large compared to the hop rate. 


D. Fast Frequency-Hopping 

There are two types of frequency-hopping: fast frequency-hopping and slow- 
frequency hopping. In fast-frequency hopping the carrier frequency changes (hops) 
more than once per data symbol and all hops within a symbol duration are combined 
in the receiver to provide the sample statistic, but in slow frequency-hopping the 
carrier frequency changes only once per data symbol. Hence, the terms of ‘slow’ 
and ‘fast’ do not mean absolutely how fast the carrier frequency changes. We define 
diversity of order L as the case where each M-ary symbol is transmitted L times over 
the channel. Thus, if /’, is the energy for each trasmission then each M-ary symbol 
requires a total energy of 


| De Od Ons (1.3) 


where £, is the symbol energy. The diversity factor is a key point of fast frequency- 
hopping communication systems, and the object of this thesis is to determine under 
what conditons fast frequency-hopping systems can provide immunity to partial-band 


interference. 


Il. ANALYSIS 


A. Determination of the Probability of Bit Error 
1. Conventional Receiver 
An analysis of the conventional FFH/MFSK receiver requires, to begin 
with, a description of the samples at the output of the M quadratic detectors be- 
fore diversity combining. The partial-band interference, when present, is assumed 
to be present in all channels. Assuming equally likely M-ary symbols, we have the 


probability of symbol error given L independent hops as 


DL 
P, = 2 ( ‘ 4'(1 — )"' P(t) (2.1) 
where P,(7) is the conditional symbol error probability given that 2 out of L hops 
are jammed. Assuming the signal is present in branch one, we have the probability 
density function of the random variable x,, representing the output of the quadratic 
detector [Ref. 3]. 


Lik + tt) fe Ga-an (aan) (2.2) 


1 
- (244/ae) = —szexp | — 
Fag (tuslae) = soa exp (TE, : 
where 2a, is the signal amplitude, 0? is noise power in hop k, and Io(e) is the 
zeroth-order modified Bessel function of the first kind. The unconditional probability 


density function is found by evaluating 


fx,,(t14) = de fx, (Tinlan) fa, (ae )dar (2.3) 
The amplitude of the signal is modeled as a Rician distributed random variable with 
a probability density function as given by (1.1). Substituting (1.1) and (2.2) into 
(2.3), we have 


eae 2a Gey 2a 
peer Coy) = x | : ( 7 


1 
~5(g?2 4902) | 2 \ 52 952 4 
roreaa (set aen) (Sea )uow) 
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The probability density function of each of the M — 1 channels with only Gaussian 


noise is 





l Umk 
IX yy (Umk) = 5g? oxP (- 55? u(rix) (2.5) 


G5. 


nt 


which is easily obtained by replacing the signal amplitude a, with 0 and replacing 7), 
with rx, Where m = 2,3,4....A7. Next, the L samples are combined to obtain the 
final sample statistic. Assuming z out of L hops are jammed, we obtain the decision 


watiagle tor branch | 
L 
Ay = a Atk 
c=) 


1 L 
ee ee (2.6) 
ke 


k>i+1 


where the superscript “) implies a jammed hop and ©) implies a non-jammed hop. 
Since each hop is assumed to be independent, we have 


Qi iy 
(2.7) 


1 2 

iD (71) ae feen(a?)| & fycola) 

where foray) is the probability density function when jamming 1s present, fyco(@3p) 
when no jamming is present, ® imphes convolution, @7 implies an 7—fold convolu- 
tion, and @(L — 2) implies an (L — 2)-fold convolution. Similarly, for the noise only 


branches, we have the random variables 


L 
Der = SS EASE 
k=1 
1 L 
= Sp eas SS x®) Wie 20, ee (2.8) 
Sa k=i+l1 


and the probability density function 
Qi @(L-1) 
EXm (Um) = Fon) @ fa (o))| 
mk mk 
For MFSk, the probability of symbol error is [Ref. 5] 
1 Hat Af—1 
PD =1= fo fry Ceski) | [ fiem(am drm] der (2.10) 
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Equation (2.10) requires fx, (x1|2) and fx,, (2m). In order to obtain fx, (z1|2), 
we first rewrite (2.7) 


fx, (wit) = fx, (zalé) ® F¥, (zalt) (22a) 


Two probability density functions are computed separately 


P(e) = [Fee] 
oF} (2.12) 


where F'y,,(s) is Laplace transform of fx,,(r1x|¢) and £7' is the inverse Laplace 


transform operation. 


LF} | = —Qajpop. (of, + os) 
sae ) a eae WE Ma eres 5 Vr (roe 5 ee 
kj kj 


] t 
a Zeles 
‘. F + 2s(of, + = ene) 


where of; represents jamming noise power plus thermal noise power. The inverse 


Laplace transform of (2.13) yields 


:—] 


2iag+z t1(o7% + 207)| ? 
1 : _ , - k 1 1 k 
ae oo | Cera] | 


tae 











11 
l 2 \/2ta?rx 

Te Tyey Jo (2.14) 

2(o4; + 207) aL; +20? 

Similarly 
L—-:1-1 
, 2(L —i)ag+a21| |21(02+207)| 2 
2 = oe ee k sot etl ees 

1 a ; 2(L —z)azr, sid 
2(o07 + 207) wb o? + 20? ( 


where of represents thermal noise power. In order to obtain fx,(z1|2), (2.14) and 


(2.15) must be convolved. In general, this cannot be done analytically. 
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For the noise branch, the probability density function 1s 


(2.16) 





} - 2 @i )  —2m @(L-1) 
€ ] | 


een oa) ars | & bet a 
k 


2G ; 


The Laplace transform of (2.16) is 


(Ole ©) ee ne 
Peale) = | Ds | ed (Qa) 


s+ po) s+ Bl”) 
where 
1 ] 
Mo, py = a 2.18 


Equation (2.17) can be expanded in a Heaviside partial fraction expansion to obtain 





on) Any Ky2 hau 
ee eae et ee ee 
s+ Be’ (s+ pf?) (s+ Bt”) 
Koy Ko KaL-i) 
* 5480 Ge a (2) (L—") Pe 
2 k ioe oO, (ee 


where the coefficients Ky;, 7 = 1,2,---,2 and Ay, , k = 1,2,---,L—12 depend on 
LE and 7. Specific coefficients for several values of L and 7 are given in Table 2.1. 


Applying the inverse Laplace transform to (2.19), we obtain 





s 2 i—1 
_ gi) : ; A432 ee 
Fxm(@m) = Px (i + Ky2tm + —S™ +--+ + i | 





2! (7-1)! 
Ko32? 
tee aa te + Ko2%m + a a 
glial | 
oe Keay] (2.20) 


We can now compute the second integral in (2.10) using (2.20) with the 


result 
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By Py 
=f, +6," 


Koy => —Niy 


L= 24 =F — 
Ky2 = (Bo)? 


Bb, - Bt 
#4 , 2 Z 
= 6 BO) = kang? — iO 


~ p46 
(3242) 
2 1 1 ‘ 
Bt — Kuby — 28 





TABLE 2.1: Coefficients of Heaviside Partial Fraction Expansion 
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ie 





| 
M- 


a4 
: (Gel ees eke 





1 
=e 
1 (1) roe) | ro -—l-r 
— Ky e777 ee 
2 Kye” GTI ny 
L-1 K 
+ i 
p=1 (P;, - 
L-i aa 


=| 
K Ly 
= Ko,e7 ——— —- —._ _--- (2.28) 
‘ La —7=1-r)i(p/y 


So the final form of P,(2) which must be evaluated numerically is 


oo Qiat tz 
At) =a - | {exp (~see Sa] 


. es 1 he | Gees }| 


2(o%; + 207) Gp, eo" 


L—-t— 
2(L —i)az+a1\ (21 (0,+207)\ ? 
QS) 4 exp 2(o7 + 207) (L —1)a? 


1 4 2(L —i)azzy 
Oe Ae ae ON 
2(a% + 207) of + 20? 








XM ior, x K a? rt a 
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2. Self-normalization Receiver 
The difference between the self-normalization receiver and the conventional 
receiver is the self-normalization factor which is the sum of all A? channel outputs 
after sampling at the outlet of the quadratic detectors. The new decision statistic is 


the combination of self-normalized quadratic detector outputs. In order to find the 


lo 


probability density function of the random variables representing the self-normalized 


quadratic detector outputs, we first define the random variable 
Vi = (Xow + Xan +-°- + Xue) —Xmk, (Mm F# 1) (2.23) 


The probability density function of the sum of independent random variables is the 
convolution of the probability density functions of those random variables. Thus, the 


probability density function of the random variable V, requires M — 2 convolutions: 


fv, (ve) 7? Fxo,(£2k) ® xa) ®+--®@ Ie (en) &) 


+ @ fxu, (Tak) (2 #m) (2.24) 


Assuming that all the probability density functions of the branches containing only 


noise are the same, we have 


fv (@r) = eee) ee = alors) (2.25) 


The Laplace transform of fx,, (xix) is 


| 


Fy,,(s) = Teel 


(2.26) 


Hence 


1 M-—-2 
Fy,(s) = era (2.27) 


Taking the inverse Laplace transform of (2.27), we get 


M-2 
_ il (Onna. VUE 
Iv, (ve) = (5 (M —3)!°P \ 903 2 
As can be seen from Fig. 1.2, the random variables that represent the self- 
normalized outputs of each of the M channels after diversity combining are obtained 


as 


L 
Us= yo Une = 1, 2 ees (2.29) 
k=] 
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Where U,;,,7 = 1,2,...,0f% & = 1,2,...,2 are the random variables that repre- 
sent the self-normalized outputs of each of the M channels for hop & of a symbol. 
These random variables are expressed in terms of the independent random variables 
Xiz,t = 1,2,...,M, &k = 1,2,...,L representing the M independent outputs of the 
M quadratic detectors for each channel of the receiver for hop & of a symbol as 


AG: 


CC — 
J Nig t+ Xan +t... + Ame 


(2.30) 


where 0 < Uj, < 1 since 0 < Xj, < 00,72 = 1,2,...,M. Using 2.29, we can express 


Pr) = U,,|2) as 


L L L 
Pr (Uy < Umi) = Pr (>> vie yD Una = Pr (>> Une = a (2.31) 


inl k= a 


Defining the random variable Z;, = Uj, — Um, we obtain 


L 
Pr(U, <Uml|i) = Pr (>: Tare 1 


a | 
= Pr(Z, <I?) (2.32) 
where 
L 
T= an (2.33) 
k=1 


is our alternative decision variable. 


From 2.30 and the definition of Z,, above, we have 


Lita ee Oe 
Nik — ee 
Nip t Noe te + Nar 
At == pea 


=k = 2.34 
Aik + Xmk + Ve (2-34) 


where —1 < Z, < 1. This random variable represents the difference between the 
power in the branch containing the signal and the power of one noise only branch 


which are detected by quadratic detectors and normalized by the self-normalization 


dey 


factor. We use the auxiliary variable method [Ref. 7] to obtain the probability density 


function of the new random variable Z,,, where we define the two auxiliary variables 
Vi, = V;, Yi = Age (2.35) 


Reforming (2.34) by replacing V, and X4, with (2.35), we have 


Y.(1 — Zin) — Zin We 


De 236 
: 1+ Zi ( ) 
and the Jacobian of the transformation jis 
(1 + 244)? 
= ——_——_—— yao 
2 ty. ( ) 


Now the joint probability density function for the random variables Z,,, W,, and 
yy, is obtained in terms of the joint probability density function of the independent 


random variables X14, Xmz, and V, as 


I 


F2in¥nW, (Z1ks Ye, Wel?) 


ee duliiieu (acne sees) fan 


| J| l + Z1k 
_ 2yk + We 
USE cap)” 
ens. ait yk + 2a; 
2(of% + 207) : 2(0? + 202) on = coal 
i —yk(1 = 216) + Z1KWe or =e 
x — exp 
2a% So ee (1 a 21k) 


1 M-2 (w,)M- =s 
‘ (a) (M3) x (- x) — 


Before evaluating the probability density function of the random variable 21;, we must 
be careful in handling the range of the variables. Note that the range of the variable 
Xmk 18 —0O < Xmxe < 00, but recall that fx, (@mx) is nonzero only for the range 
Xmk = 0. For the range of -—1 < 24, < 0, Xmz is always greater than zero; but 


for the range of 0 < Z, < 1, the range of Xm, is —co < Xm~ < oO. Consequently, 


1§ 


special care must be taken in evaluating fz,,(214). The probability density function 


of the random variable Z,; 1s obtained from 


f2,,(21k|2) ie La jaye unre | (ziks Yk, Wel? dwedy, (2.39) 


when —1 < 2; <0. Since fx,,,(¢) is nonzero only for positive values of its argument, 


then 


co pyk(1—21%)/21n 
fala) = f | jy Zanes (21s Yas Wa )ewrn dy (2.40) 


when 0 < Z,, < 1. Evaluating the integrals in (2.39) and (2.40), we have the 


probability density function of random variable Z,; conditional on 2 


a (lain —px(1 — 21) 
Iz \ 1k |?) = t a aul = —| exp (sana i ae - w) 


yee 2(1 + x) pr(1 + 1k) 
(M oa 2+ €4(1 — 21x) |: iy (ee Sb mit 


—g(z1~)u(21K) (2.41) 
where 
Bei) pay Mea avelaiea| esc cia A ( OE le zan) 
~] * a Tg a eS eee 
(1 + &,(1 — 234) )4!-2(1 + 214) 14+ & (1 — 21%) 


. l ( Pk71k )- : 
(Ce ere) ee) 
= { (Af — 2)(Af — 2)! 
+ eee: 


m=O 


: ( PZ ) 4 ( ea th ) 
een eee) nee) Le Cx zig) 


M—-2-—m 
. 221k | Pk21k | 
fr — 21,)(M —3—m)! | (1 + &(1 — 21n))(1 + &) 
Af-3-—m (= OkZ1k 
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221.(M —2—m)? 
(1 — 24)(M —2 —m-—q) 


ee el te Aue "| \I 
(1 — a) + &) 
and 
Proo= 3 
a (2.43) 


are the signal-to-noise ratios of the direct signal component and the diffuse (or faded) 
signal component, respectively. 

Having obtained the probability density function of the random variable 
Z,z, we need to consider diversity to get the overall output sample statistics. The 


random variable Z, is defined 


fe 
Zi — Aa 
k=] 
~ () (2) 
= Dy4i ar SS Lik (2.44) 
c=] k=714+1 


and the probability density function of the random variable Z, is obtained as 


Q1 @(L-i) 
(2.45) 


fz,(zi|t) = Famer) &) Foo) 


It is difficult to evaluate equation (2.45), so we leave this problem for numerical 


analysis. 


Now the probability of symbol error of the self-normalized FFH/MFSK 


receiver 1s 
Piao -Pbeaes UJ U, < U3 U oe U U, < Un|t) (2.46) 


For Af > 2, an analytic solution of (2.46) is not possible because the random variables 


U,, U2, U3,--+Unx_1 are not independent. Hence, we use the union bound to obtain 
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an upper bound on system performance. From (2.46) 
Pit) < Pr(U, < U2|2) + Pr(U, < Us|t) +---+ Pr(U, < Un|c) (2.47) 


Since we assume the sample statistics of all noise only branches are the same, then 
(2.47) becomes 


P,(i) < (M —1) Pr(Uy < Um|i) (2.48) 


where m = 2,3,4,---,M. Substituting (2.32) into (2.48), we have 


P,(i) <(M —1)Pr(Z, < 0i) (2.49) 
Since 
0 . 
Pen 0) - Ven ES (2.50) 
then from (2.49) 
O 
P,(i) < (M1) f fo, (z|i)dz1 (2.51) 


The probability of symbol error is related to the probability of bit error by [Ref. 5] 


Pi) = P,(i) x an Tina lageiit 
= esl? x aay D (2.52) 
Substituting (2.48) into (2.52), we get 
P(t) S - [- falaliaes (2.53) 


Equation (2.53) must be evaluated numerically. 


Ill. NUMERICAL ANALYSIS 


Results of the probability of error analysis for both the conventional receiver 
and the self-nomalization receiver are obtained for L = 1 to 6 hops per bit for several 
signal-to-noise ratios and direct-to-diffuse component ratios as a function of signal-to- 
interference ratio to provide an overall comparison of channel and interference effects 
on system performance for both receivers. In all cases the evaluation of either (2.22) 
or (2.41) is accomplished via a numerical integration routine. We consider L + 1 
different cases for a diversity of L hops. For example, for L = 1 the two cases to be 
considered are one non-jammed hop and one jammed hop. For L = 2 there are three 
cases: neither hop jammed, one non-jammed hop and one jammed hop, and both hops 
jammed. Additionally, the effect of the jamming ratio 7 on system performance is 
examined. We can relate the signal-to-noise ratio to the signal energy density-to-noise 


power ratio for a jammed hop as 





az + 20° - by B log, M (3.1) 
where 
ae : direct signal power component 


: diffuse signal power component 


ai +207 : total signal power 


ne : thermal noise plus jamming noise power 
ky ; average signal bit energy density 

No : thermal noise power spectral density 

Ne ; Jamming noise power spectral density 

+ : jamming ratio 


De 


B : IF bandwidth 


Te : symbol rate 
M : number of symbols in alphabet 
L : number of hops in a symbol 


For hops with no partial-band jamming 


a, +207 Ey B log, M 


oF No ie iL 





(3.2) 


where go; is the thermal noise power. Assuming the IF bandwidth (B) is equal to the 


symbol rate, we can rewrite (3.1) and (3.2) 


at + 2c? Ep log, M 


an ~ NotNj/y L 





and 


aj+20?  E,log,M 


ae No L 





For computational purposed, define 


A 
No i 
sae 
ING - 
cs PF 
202 
Og = a) 


Substituting (3.5) through (3.9) into (3.3) and (3.4), we have 


oe NYY ad 
(nt y)\(P+1)L 
ae oes nyyy W 
eq = —_ 
(n+ yv)\(p+1)L 
en 664 
0}. — 


(7 + yv) 
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Substituting (3.10), (3.11), and (3.12) into the probability density function of the 
conventional receiver and the probability density function of self-normalization re- 
ceiver and using an fast Fourier transform (FFT) routine with a sequence length of 
N=1024, numerical sequences of the probability density functions are generated for 


various signal-to-noise ratios of the jammed and non-jammed cases. 


A. Conventional Receiver 

Except for a single convolution, we have obtained analytic solution for the prob- 
ability density function of the decision statistic for the conventional FFH/MFSK de- 
tector. Unfortunately, the modified Bessel function terms make it difficult to obtain 
a complete analytic solution for the general case. Consequently, we use an FFT to 
perform the convolution numerically. Even though we have an analytic expression for 
the probability density function of the noise branch statistic (the last part of (2.22)), 
we use numerical methods to correspond with the signal statistic part which must be 
evaluated numerically. Hence, we use different equations for the numerical procedure. 


Equation (2.16) can be rewritten 


pela te in) OO) IE one) (3.13) 


where 
rm 1 Bi 
(1] Lo 7 3g3 | 
Xm(2m) = sore ks (3.14) 
20; 
z @(L—2) 
(2] es oe 
ig.) a ee 1 
Blom) = [at eo *| (3.15) 


Using the Laplace transform, we get 





(1) ae 
FU (s) = : (3.16) 
“ s+ Be 





(2) a 
faghi@y = | : =| (3.17) 


where 


1 
Oy 
1 
os as 3.19 


The inverse Laplace transform of (3.16) and (3.17) yields 





t i—l pa 50 
(1 a ! Lm 20% 3 90 
Xml tm) Ge G-1)° a) 
L-1 L-i-1 z 
ela en ee Sin ee? 3.2] 
PX (®m) (5) == See 


Now we have a different form of the probability density function of the random vari- 


able An, 





l mat peatae — 2m, 
a eed ae O.22 
(=) == / ee) 


Instead of (2.21) we substitute (3.22) into (2.22) and analyze (2.22) numerically be- 
cause we need a numerical sequence of this formula to coincide with the numerical 
sequence obtained for (2.11) The length of the basic sequence is N=1024, and prior 
to convolution it is zero-padded to create a sequence of total length N=2048 since the 
sequence is to be convolved one time. 

Numerical problems are encountered since as the value of the variable x, in- 
creases the exponential terms approach zero and the modifiied Bessel function terms 
approach infinity. Hence, the multiplication of the exponential term and the modified 
Bessel function term is invalid when the argument of the exponent is less than -745 
or the argument of the modified Bessel function is greater than 705 when 386Mat- 
lab is used. For example, in the situation described above, the value of exp(-745) is 


4.9470e-324, but the value of exp(-746) is set to 0 by the computer. The value of 
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Bessel(705),where the 386 MATLAB function for the modified Bessel function Io(e) 
is Bessel(e), is 2.2621e+304 and Bessel(706) is not defined. Even though the product 
of the two is either zero or an invalid number, the true value of a multiplication is a 
valid number. When this situation is encountered, the program evaluates the situa- 
tion and uses either zero or one (because the cumulative probability density function 
approaches one as the value of the variable increases). 

When two probability density functions are convolved by an FFT, if the valid 
range of one probability density function is much wider than the range of the other, the 
result 1s almost the same as the probability density function that has wider range. 
We must check for this situation whenever the parameters of the equation change 
while the program is running, skip the numerical convolution steps, and instead use 


the probability density function that has wider range as the convolution output. 


B. Self-Normalization Receiver 

The probability density function of self-normalization receiver (2.41) is relatively 
complicated and we must use an L-fold convolution to evaluate (2.44); an 2-fold con- 
volution for jammed hops, an (L —7)-fold convolution for nonjammed hops, and each 
result must be convolved. To implement the convolution numerically, the numerical 
sequences of the probability density function of the self-normalization receiver must 
be zero-padded to a total length of N x L before the FFT is performed in order 
to preserve the desired linear convolution of (2.44). The transformed sequences are 


multiplied point by point to produce the desired seqence of 
1 2 me 
Fz(s) = [FX(s)] x [FR(s)] (3.23) 


where fF’ is the Fourier transform of f. This result is then inverse fast Fourier trans- 
formed to yield the desired sequence for fz,(z1). The procedures for L=3 and L=5 


are similar to that for L=4 and L=6, respectively, because the length of the sequence 
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for the FFT should be a power of two, and the sequences should be zero-padded to 
iv=4096 (27°) for L=3 and P=4 and to N=6192' (2°) for L=5 and L=6. 

Equation (2.42) has many summations inside which must be carefully handled 
since numerical errors are likely when the value of the variable z;, approaches one 
When 21, approaches unity, the term 1/(1 — 2,4) approaches infinity. Therefore. we 


must group the terms properly to minimize numerical errors. Equation (2.42) is 


rewritten 


Zin(] — 216) 711 + & ace | —px(1 — 21) | 
(1+ &,(1 — zin))44-2(1 + 214) 1+ & (1 — 21x) 
é M=? (M —2)(M — 2)! | PRZ1k | 
Bane 2 ate) ce Cn — Zig) \ (le ey) 
M-3 = ea 
y aX a F Gal = =] 
221d = zr) + ef) eee Cate iat) ie a) 
(eel en) 
NaCl ie Af -—2—m Pk21k 
Vl 2 1 = 9). (ge fr eR! == 2y5)) les 7” 
i = ese Ube) spears Doerr ab) Ga) 
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in order to minimize computational errors. 


2 


IV. RESULTS 


A FFH/MFSK receiver design that demonstrates strong immunity to the partial- 
band interference is the noise normalization receiver. [’his receiver requires a measure- 
ment of the noise power per hop, but it is not easy to measure the noise power precisely 
in practical situations. This measurement is not required for the self-normalization 
receiver; hence, it is easier to implement the self-normalization receiver. As will 
be seen, partial-band jamming strongly affects the performance of the conventional 
FFH/MFSK receiver, while the self-normalization FFH/MFSK receiver performs well 
against partial-band jamming. Probability of bit error analyses for both the conven- 
tional system and the self-normalization system assume worst case partial-band jam- 
ming. The channel is modeled for a moderately strong direct-to-diffuse signal ratio 
(az/2c0? = 10) which we will refer to as Rician, and as a primarily Rayleigh faded 
channel (a%/2o0* = 0.01). For Rician fading the worst case jamming ratio typically 
increases as interference power increases. For Rayleigh fading worst case jamming 
corresponds to broadband jamming (jamming ratio y = 1) for all levels of jamming 


power. 


A. Performance of Conventional FFH/MFSK Receiver 

Graphs of the probability of bit error, including worst case, as a function of 
signal-to-jamming ratio for fixed signal-to-thermal noise ratios are obtained for 1 
through 6 hops per symbol. These results are obtained by numerically evaluating 
(2.44) for values of the jamming ratio ranging from 1.0 to 0.001 and retaining worst 
case performance for each value of E/N; . The results for a Rician fading channel 


are shown in Fig. 4.1 through Fig. 4.6. 
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Figure 4.1: The Performance of the Conventional FFH/MFSK _ nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 1, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.2: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 2, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.3: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 3, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.4: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.5: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 5, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.6: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 6, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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The signal-to-noise ratio, direct-to-diffuse ratio, and the modulation order are 
chosen to enable direct comparison with the results obtained in [Ref. 2]. We see that 
for a Rician channel that diversity is not effective in reducing the effects of worst case 
partial-band jamming. Worst case performance is poorer when diversity is used than 
for no diversity. As diversity increases, only the asymtotic value for high F,/N; ap- 
proaches the performance of the noise-nomalization receiver. The order of modulation 
has no effect with respect to providing immunity to partial-band interference. This 
is illustrated in Fig. 4.4 and Fig. 4.7 through Fig. 4.9. Generally both diversity and 
the order of modulatin are effective in improving system performance for full-band 
jamming (7 = 1) as shown Fig. 4.10 and Fig. 4.11. 

Generally both diversity and the order of modulation for a Rayleigh fading 
channel are effective in improveing system performance for full band jamming (7 = 1) 
as shown in Fig. 4.10 and Fig. 4.11. However, System performance, shown in Fig. 
4.12 through Fig. 4.16, is not improved by either diversity or modulation order when 
worst case partial-band jamming is assumed. When the signal-to-thermal noise 
ratio increases, the performance of the conventional receiver for worst case partial- 
band jamming is unchanged from that of lower signal-to-thermal noise cases except 


for the asymtotic limit of negligible Ey/N;. 


35 


Conventional FH/MFSK Receiver Performance 





E/N 213.35 dB L=4 HB 0,7/20?=10 








~ 
i 
oa 


iQ! © be pease nee 4 —— Ce ee eee : wi — 
we | © Jani ing Ratio=1.0 
ot eee hes * Jamming Ratio=8. | 
1Q° <n inten oe e ; 


ees Ng Jamming Ratio=0.01 - 
ee See ae ~Jenniing R Ratio= 0.001 






\ 
&) 
hee 
LA] 
a | 
O 1a°° \ ia \ 
| ores 
= 
i ten 
a iQ tte ae ee Pete Wh cee a glee | “octets 
) 
Solid Line : Worst base~ \ 
1Q7! 
Q 10 20 30 4Q 50 


VN ley 
b 


Figure 4.7: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 8 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.8: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 16 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.9: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 32 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.10: The Performance of the Conventional FFH/MFSK noncoher- 
ent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Jamming Ratio 


= 1, Order of Modulation = 8 and Direct-to-Diffuse Signal Power Ratio 
= 0.01 
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Figure 4.11: The Performance of the Conventional FFH/MFSK noncoher- 
ent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping Num- 


ber per Symbol = 4, Jamming Ratio = 1 and Direct-to-Diffuse Signal 
Power Ratio = 0.01 
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Figure 4.12: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 1, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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Figure 4.13: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 3, Order of Modulation = 4 and Direct-to-Diffuse 


Signal Power Ratio = 0.01 
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Figure 4.14: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 5, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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Figure 4.15: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 16 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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Figure 4.16: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 32 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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B. Performance of Self-Normalization FFH/MFSK Receiver 

All parameters used in the graphs for the self-normalization receiver are the 
same as those used for the conventional receiver to make it easier to compare with 
each other as well as with the noise-normalization receiver. As we see from Fig. 4.17 
through Fig. 4.25, for a Rician fading channel, increasing diversity and increasing 
modulation order are both effective in minimizing the effects of worst case partial- 
band interference. 

Overall, the performance of the self-normalization receiver is almost the same 
as that of the noise-normalization receiver. There exists a small difference in the 
performance of the two receivers. This difference is expected due to the union bound 
used to obtain the FFH/MFSK self-normalized receiver performance. The union 
bound has the most pronounced effect on system performance when E£;/N; is small 
and the modulation order is high. See, for example, Fig. 4.26 where the probability 
of bit error for a bit energy-to-jamming noise density ratio of 0 dB with a jamming 
ratio 1.0 is shown greater than unity. This is not valid for a probability. However, 
the error performance at low E/N; is not important, and the error due to the union 
bound compared to the exact value at high F,/N; is not significant. 

The performance of the self-normalization receiver for a Rayleigh fading channel 
in worst case partial-band jamming is the same as for full band jamming (y = 1) 
regardless of the parameters as shown in Fig. 4.27 through Fig. 4.29. Generally for 
a Rayleigh fading channel, diversity is effective in improving the performance of the 
self-normalization receiver. This is shown in Fig. 4.30. The effect of increasing the 
modulation order is less significant as shown in Fig. 4.26. As mentioned previously, 


the union bound makes the simulated performance worse than the true performance. 
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Figure 4.17: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 1, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.18: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 2, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.19: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 3, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.20: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.21: The Performance of the Self-Normalization FFH/MESK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 5, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.22: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol! = 6, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.23: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 8 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.24: The Performance of the Self-Normalization FFH/MEFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 16 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.25: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 


Number per Symbol = 4, Order of Modulation = 32 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.26: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping 


Number per Symbol = 4, Jamming Ratio = 1 and Direct-to-Diffuse Signal 
Power Ratio = 0.01 
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Figure 4.27: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0 dB, Hopping 


Number = 1, Order of Modulation = 4 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 
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Figure 4.28: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping 


Number = 4, Order of Modulation = 4 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 
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Figure 4.29: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping 


Number = 6, Order of Modulation = 32 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 
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Figure 4.30: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Jamming 


Ratio = 1, Order of Modulation = 8 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 
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V. CONCLUSION 


The probability of bit error performance for a conventional noncoherent FFH/MFSk 
receiver has been obtained for Rayleigh- and Rician-faded channel with partial- 
band jamming. The effect of fading is detrimental for both fading models, but fast 
frequency-hopping in general provides a means to overcome, at least partially, fading 
effects. However, the conventional receiver is severely affected by worst case partial- 
band jamming. It does not show any improvement in worst case performance with 
high diversity or high order of modulation. The probability of bit error of the conven- 
tional receiver in worst case partial-band jamming is always between 107? and 107° 
at E,/N; of 20 dB regardless of E,/No, the order of modulation, diversity factor, or 
direct-to diffuse power ratio. The probability of bit error of the self-normalization re- 
ceiver in worst case partial-band jamming is the same as that of conventional receiver 
for slow frequency-hopping (L=1). However, for fast frequency-hopping, the perfor- 
mance becomes better and better with increasing diversity factors: the probability of 
Prmecroris) ~ 109" at L=2- 2 « 10m at G=3-10-- at L=4 and 8 107° at L=5 and 
L=6. When the diversity factor is 6, the worst case performance is almost same as 
that of full band jamming case (jamming ratio 1.0). At this point the jammer does 
not need any information about the communicator to optimize the jamming effect by 
managing the jamming power. 

The self-normalization receiver with diversity provides very good immunity to 
worst case partial-band jamming and has a performance similar to that of the noise- 
normalization receiver. The order of modulation also improves the overall perfor- 
mance of the self-normalization receiver, but partial-band jamming affects higher 


orders of modulation more than the lower ones. Thus, if the communicator uses a 


6] 


higher order of modulation and wants to suppress all partial-band jamming effects, 


he should use a higher hopping rate. 
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APPENDIX A 


This is the Matlab program to produce the probability of bit % 


error of 


% 
% 


ee =— Ps OO 


oe 


dP de dO oe dP dO de oe de 


Conventional FH/MFSK quadratic receiver. 


Lee,Kang Yeun (Major ROKAF) 

clear 

3 signal-to-noise ratios to be implemented (13.35 db, 16 db and 
20 db) 

SNR=[21.62718524D0 39.81071706D0 100.0D0]}; 

2 direct-to-diffuse signal power ratios (0.01 and 10) 
DDR=[0.01D0 10.0D0]; 

4 jamming ratio (1, 0.1, 0.01 and 0.001) 


JRATIO=[1.0D0 0.1D0 0.01D0 0.001D0}; 
data=zeros(8,51); 


set the thermal plus jamming noise power 1 
rhokl=1; 

basic sequence length is 1024 

lim=1024; 

n=2*lim; 


zpad=zeros(1,lim); 

opad=ones(1,1lim); 

hopping number (1 to 6) 

for L=1:6, 

for j1=1:2, psi=DDR(j1); 

for }2=1:3, nu=SNR(3j2); 

for }33=1:4, gama=JRATIO(j3); 

the order of modulation (M-ary : 4, 8, 16, 32-ary) 
for }4=2:5; M=2%34; 

for 35=0:2:50, i1ota=10.0D0%*(0.1D0*j5) ; 

diffuse signal power ( rho ) 

rho=nu*gama*iota/ (nu+gama*iota) /(pSit1) *}4/L; 
direct signal power ( alpha ) 

alpha=rho*psi; 

thermal noise power ( rhok2 ) 

rhok2=iota*gama/ (iota*gama+nu) ; 

set the sequence limit roughly based on the singal-to-noise ratio 
(nu) 

limx=80*nu ; 

Givide the range by 100 to check the proper limit 
delta=limx/100; 

select almost the middle value of the increment to compute the 
sequence 

of the probability density function. 
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0 dO dO 


0 dO dO 


x=0.49*delta: delta: limx-0.51*delta; 
g=exp (-(x+2*L*alpha) /2/ (rhokl1+rho) )/2/ (rhokl+rho). 

* x/2/L/alpha).%... 

((L-1)/2) .*abs(j* (L-1) *bessel (L-1,j*sqrt(2*L*alpha*x) 

/ (rhoki1+rho) ))*delta; 
check the sequence whether it has invalid sequence or not 
fi=find(isnan(g) ); 
g(f1)=zeros(1,length(fi)); 
if the sequence has only Os, the limit is too wide. Thus set the 
first 
increment as the new limit 
sez all(g==0)==1, 
gomin=1; 
if the sequence is valid, check the point around 1.0e-6 to select 
the 
new limit 
else, [xx,gmax]=max(g) ; 
g(1:gmax)=ones(1,gmax) *xx; 
(xx, gmin]=min((g-1.0e-6) .*2); 
end 
limx=gmin*delta; 
if the new limit is first or second value of the increment 
redo above procedure 
while gmin <= 2, 
delta=limx/100; x=0.5*delta:delta:limx-0.5*delta; 
g=exp (- (x+2*L*alpha)/2/(rhokl+rho))/2/ ... 

(rhokl+rho) .*(x/2/L/alpha).* ... 

((L-1)/2) .*abs(j* (L-1) *bessel (L-1,j}*sqrt(2*L*alpha*x) ... 

/ (rhnoki1+rho) ))*delta; 
fi=find(isnan(g) );? 
g(f1)=zeros(1,length(fi)); 
if all (g==0) ==1, 
gmin=1; 
else, [xXx,gmax]=max(g); 
g(1:gmax) =ones(1,9gmax) *xx; 
[xx,9gmin]=min((g-1.0e-6) .*2); 
end 
limx=gmin*delta; 
end 
after setting the limit, whole range is divided by 1024 
to produce the sequence 
delta=limx/lim;x=0.49*delta:delta:1limx-0.51*delta; 
if jamming ratio is 1, the computation is simpler than others 
if 43 == 1, 
gs=- (x+2*L*alpha) /2/ (rhok1+rho) -log(2* (rhok1+rho) ) 
+log(x/2/L/alpha)* ... 
(L-1) /2+log (abs (j* (L-1) *bessel(L-1,j*sqrt(2*L*alpha*x) ... 
/ (rhok1+rho) )));3 
gs=exp(gs) *delta; 
fi=find(isnan(gs) ); 
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gs(fi)=zeros(1,length(fi)); 

if all(gs==0)==1, 

gs=zeros(1,lim); 

gs(1)=1; 

else,gs=gs/sum(gs); 

end 
ps=log(1/2/rhok1) *L4+log(x) *(L-1)-log(fact(L-1)) - x/2/rhokl; 
ps=exp (ps) *delta; 

fi=find(isnan(ps)); 

ps (fi)=zeros(1,length(fi)); 

if all(ps==0)==1, 

ps=ones(1,1lim); 

else,ps=ps/sum(ps); 

ps=cCumsum(ps) ; 

ps=ps/max (ps) } 

end 

data (count, j5+1)=M/2/ (M-1) * (1l-sum(gs.*ps.*(M-1))); 
else 

if jamming ratio is not 1 , the hopping number has to be 
considered 

For 1 = 0: L 

if no hop is jammed 

if 1 == 0, 

delta=limx/100; 

x=0.5*delta: delta: limx-0.5*delta; 

g=- (x+2*L*alpha) /2/ (rhok2+rho) -log(2* (rhok2+rho) ) 

+log(x/2/L/alpha) * ‘ 

(L-1) /2+log (abs (j* (L-1) *bessel (L-1,j}*sqrt(2*L*alpha*x) 

/(rhok2+rho)))); 

g=exp(g) *delta; 

fi=find(isnan(g) ); 

g(fi)=zeros(1,length(fi)); 

if all (g==0)==1, 

gmin=1; 

else, [xx,gmax]=max(g) 7;g(1:gmax)=ones(1,gmax) *xx; 
(xx,gmin]=min((g-1.0e-6) .*2); 

end 

vlimx=gmin*delta; 

check the sequence limit again 

while gmin <= 2, 

delta=vlimx/100; 
x=0.5*delta:delta:vlimx-0.5*delta; 

g=- (x+2*L*alpha) /2/ (rhok2+rho) -log(2* (rhok2+rho) ) 
+log(x/2/L/alpha) * ‘ 

(L-1) /2+log (abs (j* (L-1) *bessel (L-1,j)*sqrt(2*L*alpha*x) 
/(rhok2+rho)))); 

g=exp(g) *delta; 

fi=find(isnan(g) ); 

g(f1)=zeros(1,length(fi)); 


65 


1f all (g==0)==1, 

gmin=1; 

else, 

[xx,gmax]=max(g) ;g(1:gmax) =ones (1, gmax) *xx; 

(xx, gmin]=min((g-1.0e-6) .*2); 

end 

vlimx=gmin*delta; 

end 

delta=vlimx/lim; 

x=0.49*delta:delta:vlimx-0.51*delta; 

gs=- (x+2*L*alpha) /2/ (rhok2+rho) -log(2*(rhok2+rho)) ... 
+log(x/2/L/alpha)* ... 

(L-1) /2+log(abs(j* (L-1) *bessel (L-1,j*sqrt(2*L*alpha*x) 

/ (rhok2+rho) )));3 

gs=exp (gs) *delta; 

fi=find(isnan(gs) ); 

gs(fi)=zeros(1,length(fi)); 

1f all(gs==0)==1, 

gs=zeros(1,lim); 


end 
ps=-log(2*rhok2) *L+log(x) *(L-1)-log(fact(L-1)) - x/2/rhok2; 
ps=exp (ps) *delta; 
fi=find(isnan(ps) ); 
ps(fi)=zeros(1,length(fi)); 
if all(ps==0) ==1, 
ps=ones(1,1lim) ; 

else, ps=ps/sum(ps) ; 
ps=ps/sum(ps) ; 
ps=cumsum(ps) ; 
PpSs=ps/max (ps) ; 


end 
if all hops are jammed 
elseif i == L, 


delta=limx/100; 

x=0.5*delta:delta: limx-0.5*delta; 

g=- (x+2*1i*alpha) /2/ (rhok1+rho) -log(2*(rhokli+rho)) ... 
+log(x/2/i/alpha)* ... 

(1-1) /2+log(abs(j%* (i-1) *bessel(i-1,j*sqrt(2*i*alpha*x) ... 
/(rhoki+rho))));3 

g=exp(g) *delta; 

fi=find(isnan(g)); 

g(f1)=zeros(1,length(fi)); 

if all(g==0)==1, 

gmin=1; 

else, 

[xx,gmax]=max(g) ;g(1:gmax) =ones(1,gmax) *xx; 
[xx,gmin]=min((g-1.0e-6) .*2); 
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end 

ulimx=gmin*delta; 

check the sequence limit again 

while gmin <=2, 

delta=ulimx/100; 

x=0.5*delta:delta:ulimx-0.5*delta; 
g=-(x+2*i*talpha) /2/ (rhokl+rho) -log(2*(rhokl+rho)) ... 
+log(x/2/i/alpha)* ... 
(i-1)/2+log(abs(j*(i-1) *bessel (i-1,j*sqrt(2*1i*alpha*x) 
/ (rhokl+rho)))); 

g=exp(g) *delta; 

fi=find(isnan(g)); 

g(fi)=zeros(1,length(fi)); 

if all(g==0)==1, 

gomin=1; 

else, 

(xx,gmax]=max(g) ; 

g(1:gmax) =ones(1,gmax) *xx; 

[xx, gminJ=min((g-1.0e-6) .*2); 

end 


ulimx=gmin*delta; 


end 

delta=ulimx/lim; 

x=0.49*delta:delta:ulimx-0.51*delta; 

gs=- (x+2*i*alpha) /2/ (rhokl1+rho) -log(2* (rhok1+rho) ) 
+log(x/2/i/alpha)* ... 

(i-1) /2+log (abs (j*(i-1) *bessel (i-1,j*sqrt(2*1*alpha*x) 


/(rhoki+rho)))); 

gs=exp(gs) *delta; 
fi=find(isnan(gs)); 
gs(fi)=zeros(1,length(fi)); 
if all(gs==0)==1, 
gs=zeros(1,lim); 

gs(1)=1; 

else, 

gs=gs/sum(gs) ; 

end 
ps=-log(2*rhok1) *itlog(x) *(i-1)-log(fact(i-1)) - x/2/rhokl; 
ps=exp (ps) *delta; 

if all(ps==0) ==1, 

ps= ones(1,lim); 

else, 

fi= find(isnan(ps) ); 

ps (f1)=zeros(1,length(fi)); 
ps = cumsum(ps); 

pS = psS/max(ps) } 

end 
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else 

if i hops are jammed and L-i hops are not jammed 
check the sequence limit again 

rhok=rhok1; 

delta=limx/100; 

x=0.5*delta:delta:limx-0.5%*delta; 

g=- (x+2*1i*alpha) /2/ (rhok+rho) -log(2*(rhok+rho)) ... 
+tlog(x/2/i/alpha)* ... 
(i-1)/2+log(abs(j* (i-1) *bessel (i-1,)*sqrt(2*i*alpha*x) 
/(rhok+rho))))?; 

g=exp(g) *delta; 

fi=find(isnan(g) ); 

g(fi)=zeros(1,length(fi)); 

if all (g==0)==1, 

gmin=1; 

else, 

[xx,gmax]=max(9g) ; 

g (1: gmax) =ones (1,9gmax) *xx; 
(xx,gmin]=min((g-1.0e-6) .*2); 

end 


wlimx=gmin*delta; 

while gmin <=2, 

delta=wlimx/100; 

x=0.5*delta:delta:wlimx-0.5*delta; 

g=- (x+2*i*alpha) /2/ (rhok+rho) -log(2*(rhok+rho)) ... 

+log(x/2/i/alpha)* ... 

(1-1) /2+log(abs(j*(i-1) *bessel (i-1,j*sqrt(2*i*alpha*x) ... 

/ (rhok+rho) ))); 

g=exp(g) *delta; 

fi=find(isnan(g) ); 

g(f1)=zeros(1,length(fi)); 

if all(g==0)==1, 

gmin=1; 

else, 

[xx,gmax]=max(g) ; 

g(1:gmax) =ones(1,9gmax) *xx; 
(xx,gmin]=min((g-1.0e-6) .*2); 

end 

wlimx=gmin*delta; 

end 

rhok=rhok2; 

delta=1limx/100; 

x=0.5*delta:delta:limx-0.5*delta; 

g=- (x+2*(L-1) *alpha) /2/ (rhok+rho) -log(2*(rhok+rho) ) 
+log(x/2/(L-1) /alpha) * 

(L-1-1) /2+log(abs(j* (L-i-1) *bessel (L-i-1,j*sqrt(2*(L-1) ... 
*alpha*x)/ ... 

(rhok+rho) ))); 
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g=exp(g) *delta; 
fi=find(isnan(g) ); 
g(fi)=zeros(1,length(fi)); 
if all(g==0)==1, 

gmin=1; 

else, 

[xx,gmax]=max(g) ; 
g(1:gmax)=ones (1, gmax) *xx; 
(xx,gmin]=min((g-1.0e-6) .*2); 
end 

zlimx=gmin*delta; 


check the sequence limit again 

while gmin <=2, 

delta=zlimx/100; 

x=0.5*delta:delta:zlimx-0.5*delta; 

g=- (x+2* (L-1) *alpha) /2/ (rhok+rho) -log(2* (rhok+rho) ) 
+log(x/2/(L-i1)/alpha)* ... 

(L-i-1) /2+1log (abs (j* (L-i-1) *bessel (L-i-1,j*sqrt(2*(L-1i) 
*alpha*x)/(rhok+rho)))); 

g=exp(g) *delta; 

fi=find(isnan(g) ); 

g(f1)=zeros(1,length(fi)); 

if all(g==0)==1, 

gmin=1; 

else, 

(xx,gmax])=max(g); 

g(1:gqmax)=ones (1,gmax) *xx; 

(xx, gmin]=min((g-1.0e-6) .*2); 

end 

zlimx=gmin*delta; 

end 

produce the sequency of jammed hop of signal branch 
delta=wlimx/lim; 
x=0.49*delta:delta:wlimx-0.51*delta; 
gl=-(x+2*1i*alpha) /2/ (rhok1+rho) -l1o0g(2* (rhokl+rho) ) 
+tlog(x/2/i/alpha)* ... 

(1-1) /2+log(abs(j*(i-1) *bessel (i-1,j*sqrt(2*i*alpha*x) 
/(rhokl+rho)))); 

gl=exp(gl) *delta; 

fi=find(isnan(gl) ); 

g1(f1)=zeros(1,length(fi)); 

if all(gl==0)==1, 

gl=zeros(1,lim); 

g1(1)=1; 

else, 

gl=gl1/sum(gl) ; 

end 

produce the sequence of nonjammed hop of signal branch 
if wlimx >= 100*zlimx 
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gs=(gl zpad]; 

else 

g2=-(x+2* (L-1) *alpha) /2/ (rhok2+rho) -log(2* (rhok2+rho) ) 
+log(x/2/(L-1) ... 

/alpha) *(L-i-1) /2+log(abs(j*(L-i-1) *bessel(L-i-1, ... 
jJ*sqrt(2*(L-i) ... 

*alpha*x) /(rhok2+rho) ))); 

g2=exp(g2) *delta; 

fi=find(isnan(g2)); 

g2(fi)=zeros(1,length(fi)); 

if all(g2==0)==1, 

g2=zeros(1,lim); 

g2(1)=1; 

else, 

g2=9g2/sum(g2) ; 

end 


produce the sequence of signal branch 
gs = abs(ifft(fft(gl,n).*fft(g2,n),n)); 
gs = gs/sum(gs) ; 

end 


produce the sequence of jammed hop of noise branch 
pl=-log (2*rhok1) *i+log(x) *(i-1) -log(fact(i-1) )-x/2/rhok1; 
pl=exp(pl1) *delta; 

fi=find(isnan(pl)); 

p1(fi)=zeros(1,length(fi)); 

if all(pl==0)==1, 

pl=zeros(1,lim); 

p1(1)=1; 

else, 

Pl=pl/sum(p1) ; 

end 

produce the sequence of nonjammed hop of noise branch 
p2=-1log (2*rhok2) * (L-1) +log(x) * (L-i-1) -log(fact (L-1i-1) ) -x/2/rhok2; 
p2=exp(p2) *delta; 

fi=find(isnan(p2)); 

p2(f1)=zeros(1,length(fi)); 

if all(p2==0)==1, 

p2=zeros(1,lim); 


p2(1)=1; 

else, 
p2=p2/sum(p2) ; 
end 


produce the sequence of noise branch 


ps = abs(ifft(fft(pl,n).*fft(p2,n),n)); 
ps = cumsum(ps) ; 

ps = ps/max(ps) ; 

end 
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% compute the probability of bit error using sequence 
tdata (j5+1)=tdata(j5+1)+bino(L,i) *gama*i* (1-gama) * (L-1i) 
*M/2* (l-sum(gs.*ps.*(M-1))); 
end, end, end 
data=[data;td]; 
save datafile data; 


end,end,end,end,end 


de de de GP AP AO AO 0 dO dP 


APPENDIX B 


This is the program to produce the probability of bit error of 
Self-Normalization FH/MFSK Receiver 
Lee, Kang Yeun (Major, ROKAF) 


three special routines 

fact : factorial routine 

bino : binomial coefficient routine 

ereval : error probability computation routine using fft 


3 signal-to-noise ratio(13.35db, 16db, 20db) 
SNR=[21.62718524D0 39.81071706D0 100.0D0]}; 

3 direct-to-diffuse signal power ratio (0.01, 1, 10) 
DDR=[0.01D0 1.0D0O 10.0D0}; 

4 jamming ratio 

JRATIOI=[1.0D0 0.2D0 0.1D0 0.05D0 0.01D0]; 
JRATIO2=[1.0D0 0.2D0 0.05D0 0.01D0 0.001D0}; 


Bandwidth-to-Symbol Rate Ratio (BW/Rs : 1 2 ) 
mu=1.0D0; 

M-ary 

M=16.0D0; 

Multichannel factor w=logM/log2 


omega=log(M) /log(2.0D0) ; 


Bit SNR 13.35 db 
nu=SNR(1) ; 


count=1; 


Hopping number (1 to 6) 
for J1=1:6 
L=J1; 


Direct-to-Diffuse signal ratio ( 0.01 1 10 ) 
for J3=1:3 

psi=DDR(J3) ; 

Direct SNR of NOJAMMING CASE 
rho=omega*nu*psi/mu/ (psit+1.0D0O) /L; 

Diffuse SNR of NOJAMMING CASE 

Xi=omega*nu/mu/ (psi+1.0ODO) /L; 


produce the sequence of nojamming case 
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z=-0.999D0:0.002D0:0.999D0; 

NOJAM=exp (rho*(z-1.0D0)./(2.0D0+xi*(1.0D0-z))). ... 
*(1.0D0+z).*(M-2.0D0O)... 

./(2.0D0+xi*(1.0D0-z) ).*((2.0D0*(1.0D0+xi))./ ... 

tee er lec ©. 0) se St 2 OD O> (DP iv ert 
.*(1.0D0+rho*(1.0D0+z) ./((2.0D0+xi*(1.0D0-z)) ... 
*(1.0D0+xi) ) )+M-2.0D0); 


z=0.001D0:0.002D0:0.999D0; 
B=rho*z./((1.0D0+xi*(1.0D0-z) )*(1.0D0+x1) ); 
E=exp (rho*(z-1.0D0) ./(1.0D0+x1i*(1.0D0-z))); 


TT=zeros (M-1,500); 

for N=0:M-2 

TT (N+1,:)=fact (M-2.0D0)/fact (M-2.0DO0-N) /fact(N)*2*B.“N; 
end 

Tl=sum(TT) .*((1.0D0-z) *(1.0D0+x1)).*(M-2).*z. ... 
*(M-2.0D0)./(1.0D0+z)... 

-/(1.0D0+x1*(1.0D0-z) ).*(M-1).*E; 


TT=zeros (M-2,500); 

for N=0:M-3 

TT1=zeros(M-1-N,500) ; 

for IQ=0:M-2-N 

Primi Orl, :)—fLact (M-2.0D0-N) /fact (M-2. ODO-N-10) / fact (10) *2*B.*10; 

end 

TEMP=2.0D0*((1.0D0-z) *(1.0D0+xi) ) .*(M-3-N) .*(1.O0D0+x1) .*z2.%*2 ... 
~*(1.0D0+z) .*(N-1) .*(M-2.O0DO-N) ; 

TEMP=TEMP./(1.0D0+x1i*(1.0D0-z2) ).*(M-1.ODO-N) .*E.*sum(TT1) ; 

TT1=zeros (M-2-N,500) ; 

for IQ0=0:M-3-N 

TT1(IQ+1,:)=fact(M-3.O0DO-N) /fact (M-3.ODO-N-IQ) /fact (IQ) *2*B.*I1Q; 

end 

TEMP1=((1.0D0-z) *(1.0D0+x1) ).*(M-3-N) .*z2.*(1.0D0+2) ... 
~“N.*(M-2.0D0) ./(1.0D0+x1* (1.0D0-z2)) .* (M-2-N) .*E; 

Liekos—=— 0) 

TT (N+1,:)=TEMP+TEMP1.*TT1; 

else 

TT (N+1,:)=TEMP+TEMP1.*sum(TT1) ; 

end 

end 

T1=T1+sum(TT) ; 


% sequence of nojamming case 
NOJAM (501:1000)=NOJAM(501:1000)-T1; 


Ps=zeros(5,41); 


% jamming ratio 
for J4=1:5 


if J3== 
gamma=JRATIO2 (J4) ; 
else 
gamma=JRATIO1(J4); 
end 


% Signal-to-Jamming noise ratio ( 0 db - 40 db ) 


for J5=0:40 
1ota=10.0D0% (0.1D0*J5) ; 


% Direct SNR of JAMMING case 
rho=omega*nu*iota*psi*gamma/mu/ (psit+1.0D0) / (lota*gamma+nu) /L; 

% Diffuse SNR of JAMMING case 
x1l=omega*nu*iota*gamma/mu/ (psi+1.0D0) / (iota*gamma+nu) /L; 


% produce the sequence when the signal is jammed by noise jamming 


Z=-0.999D0:0.002D0:0.999D0; 
JAM=exp (rho*(z-1.0D0) ./(2.0D0+x1i*(1.0D0-z))). ... 
*(1.0D0+z).*(M-2.0D0)... 

-/(2.0D0+xi*(1.0D0-z)). ... 
*( (2.0D0*(1.0D0+xi))./(xi*(1.0D0-z)+2.0D0)... 
-*(1.0D0+rho*(1.0D0+z) ./((2.0D0+x1*(1.0D0-z)) ... 
*(1.0D0+xi) ))+M-2.0D0) ; 


Z=0.001D0:0.002D0:0.999D0; 
B=rho*z./((1.0D0+xi*(1.0D0-z) )*(1.0D0+x1) ); 

E=exp (rho*(z-1.0D0) ./(1.0D0+xi*(1.0D0-z))); 

TT=zeros (M-1,500) ; 

for N=0:M-2 

TT (N+1,:)=fact (M-2.0D0) /fact (M-2.O0DO-N) /fact (N) *2*B.“N; 
end 

T1l=sum(TT) .*((1.0D0-z) *(1.0D0+x1i)).*(M-2). ... 
*Z2.*(M-2.0D0)./(1.0D0+zZ)... 
-/(1.0D0+x1i*(1.0D0-z) ).*(M-1) .*E; 


TT=zeros (M-2,500) ; 

for N=0:M-3 

TT1=zeros (M-1-N, 500) ; 

for IQ=0:M-2-N 

TT1(IQ+1,:)=fact (M-2.O0D0-N) /fact (M-2.0DO-N-IQ) /fact(1IQ) *2*B.*IQ; 

end 

TEMP=2.0D0*((1.0D0—-z) *(1.0D0+xi) ) .*(M-3-N) .*(1.0D0+x1) .*z.*%2 ... 
~-*(1.0D0+Z) .*(N-1) .* (M-2.0DO=N) ; 

TEMP=TEMP./(1.0D0+x1i*(1.0D0-z) ).*(M-1.O0DO-N) .*E.*sum(TT1) ; 

TT1=zeros (M-2-N,500) ; 

for I0=0:M-3-N 

TT1(IQ+1,:)=fact (M-3.0DO-N) /fact (M-3.0DO-N-IQ) /fact (IQ) *2*B.*1IQ; 
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end 

TEMP1=((1.0D0-z) *(1.0D0+xi) ) .* (M-3-N) .*Z.*(1.0D0+2z) 
eon. *(M—-2.0D0) ./ (1.0D0+x1i* (1.0D0-Z) ) .*(M-2-=N) .*E; 

if IQ == 0, 

TT (N+1,:)=TEMP+TEMP1.*TT1; 

else 

TT (N+1, :)=TEMP+TEMP1.*sum(TT1) ; 

end 

end 

T1l=T1l+sum(TT) ; 


% the sequence of jamming case 
JAM (501:1000)=JAM(501:1000) -T1; 


% 'bino' and ‘ereval' are the functions for binomial and error 
evaluation. 


% 

% Ps(i) <= (M-1) * Pr(z<0O) 

% Pb = (M/2)/(M-1)*Ps 

% Pb <= (M/2) * Pr( z< oO) 

% probability of symbol error 

for 1=0:L 

Ps (J4,J35+1)=Ps (J4,J35+1)+bino(L,i) *gamma‘*i* (1-gamma) * (L-1) 
*ereval (NOJAM,L-1,JAM,1); 

end 


% probability of bit error 
Ps (J4,J35+1)=(M/2.0D0) *Ps(J4,35+1); 


end 
end 


% Save the Results and Plot 
save data Ps 


end 
end 
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TUE 


ay. 
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